mm duq 

HUMAN NEUROSCIENCE 



ORIGINAL RESEARCH ARTICLE 

published: 23 July 2014 
doi: 10.3389/fnhum.2014. 00529 




LORETA EEG phase reset of the default mode network 

Robert W. Thatcher*, Duane M. North and Carl J. Biver 

EEG and Neurolmaging Laboratory, Applied Neuroscience Research Institute, Seminole, FL, USA 



Edited by: 

Jing Xiang, Cincinnati Children's 
Hospital Medical Center, USA 

Reviewed by: 

Arun Bokde, Trinity College Dublin, 
Ireland 

Xiaoshan Wang, Nanjing Brain 
Hospital, China 

'Correspondence: 

Robert W. Thatcher, EEG and 
Neurolmaging Laboratory, Applied 
Neuroscience Research Institute, 
7985 113th Street, Suite 210, 
Seminole, FL 33772, USA 
e-mail: rwthatcher2@yahoo.com 



Objectives: The purpose of this study was to explore phase reset of 3-dimensional current 
sources in Brodmann areas located in the human default mode network (DMN) using Low 
Resolution Electromagnetic Tomography (LORETA) of the human electroencephalogram 
(EEG). 

Methods: The EEG was recorded from 19 scalp locations from 70 healthy normal subjects 
ranging in age from 13 to 20 years. A time point by time point computation of LORETA 
current sources were computed for 14 Brodmann areas comprising the DMN in the delta 
frequency band. The Hilbert transform of the LORETA time series was used to compute 
the instantaneous phase differences between all pairs of Brodmann areas. Phase shift and 
lock durations were calculated based on the 1st and 2nd derivatives of the time series of 
phase differences. 

Results: Phase shift duration exhibited three discrete modes at approximately: (1) 25 ms, 
(2) 50 ms, and (3) 65 ms. Phase lock duration present primarily at: (1) 300-350 ms and 
(2) 350-450 ms. Phase shift and lock durations were inversely related and exhibited an 
exponential change with distance between Brodmann areas. 

Conclusions: The results are explained by local neural packing density of network hubs 
and an exponential decrease in connections with distance from a hub. The results are 
consistent with a discrete temporal model of brain function where anatomical hubs behave 
like a "shutter" that opens and closes at specific durations as nodes of a network giving 
rise to temporarily phase locked clusters of neurons for specific durations. 
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INTRODUCTION 

When one is at rest and not engaged in a task and absorbed in 
a ruminating self-narrative about the past and future then it is 
during these reflective moments that the default mode network 
(DMN) is activated and the attention network is anti-correlated 
or reciprocally deactivated (Raichle et al., 2001; Raichle, 2010). 
The insula appears to function as a switch that is correlated with 
phase shifting of the attention and default networks activation 
vs. suppression (Bressler and Menon, 2010). Petersen and Posner 
(2012) review the functional MRI (fMRI) studies of the atten- 
tion network and the DMN in attention deficit disorders that are 
characterized by the intrusion of the self-narrative in academic 
situations resulting in poor grades. The reciprocal relationship 
between the DMN related to an ongoing internal self-narrative 
and the attention network focused on the external world is an 
important dynamic, however, fMRI has a limited temporal res- 
olution and is unable to resolve millisecond periods of phase 
lock and phase shift of neurons located in network nodes and 
functional connections that comprise the DMN. 

The EEG has slightly less spatial resolution than fMRI, but ade- 
quate spatial resolution to measure the average current density 
of Brodmann areas at 2 cm 3 to about 3 cm 3 volumes in the mil- 
lisecond time domain (Pascual-Marqui, 1999; Vitacco et al., 2002; 
Mulert et al., 2004; Grech et al., 2011). 



There are changes in the synaptic synchrony of millions of 
neurons connected at varying time delays and frequencies in 
the DMN. The "DMN" is constituted primarily by the cingu- 
late gyrus, hippocampus, medial frontal lobes, temporal lobes 
and parietal lobes with approximately five times the number of 
synaptic connections than any other cortical network (Buckner 
et al., 2008; Hagmann et al., 2008). Activation of the DMN signif- 
icantly increases demand on blood glucose and oxygen as well as 
changes in the synchrony of synaptic potentials on the dendrites 
and cell bodies of cortical pyramidal neurons as measured in the 
human EEG using 3-dimensional electrical neuroimaging meth- 
ods (Pascual-Marqui et al., 1994; Pascual-Marqui, 1999; Michel 
et al., 2009) also referred to as EEG Tomography (tEEG) (Cannon 
et al., 2009; Thatcher, 2011; Thatcher et al., 2011) or Brain 
Electromagnetic Tomography (BET) (Valdes-Sosa et al., 1992; 
Bosch-Bayard et al., 2001; Hernandez-Gonzalez et al., 2011). 

Hughes and Crunelli (2007) and Buzsaki (2006) review how 
action potentials occur when neurons are in-phase with respect to 
the local field potentials (LFPs) and how action potentials are sup- 
pressed when neurons are shifted anti-phase with respect to the 
LFP where phase shift is a high speed switch for large collections 
of neurons to functionally synchronize with sub-sets of neurons 
in different network nodes in the millisecond time domain. The 
human EEG is the summation of LFPs arising from pyramidal 
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neuron synapses. The process of phase shift and phase lock of 
the EEG between different scalp locations has been shown to be 
correlated with a range of functional conditions. For example, 
measures of EEG phase reset (PR) have been correlated to vari- 
ous frequency bands during cognitive tasks (Tesche and Karhu, 
2000; Kirschfeld, 2005; Kahana, 2006), working memory (John, 
1968; Damasio, 1989; Tallon-Baudry et al., 2001; Rizzuto et al, 
2003), sensory-motor interactions (Vaadia et al, 1995; Roelfsema 
et al., 1997), hippocampal long-term potentiation (McCartney 
et al., 2004), brain development (Thatcher et al., 2008b), autism 
(Thatcher et al, 2009b), and consciousness (Varela et al., 2001; 
John, 2002, 2005; Cosmelli et al, 2004). 

These studies indicate that measures of phase shift duration 
and phase lock duration between groups of neurons measured 
at the scalp EEG must also be capable of being measured at the 
level of the 3-dimensional sources of the EEG using inverse meth- 
ods. The initiation of in-phase to anti-phase dynamics of neurons 
is modeled by long distant excitatory inputs on dendrites where 
phase shift duration is dependent on the density of neurons in 
local loops and phase lock duration is determined by the reaction 
to the short duration excitatory inputs producing long duration 
inhibitory synaptic potentials (van Drongelen et al, 2004; Ko and 
Ermentrout, 2007; Tiesinga and Sejnowski, 2010; Li and Zhou, 
2011). These studies suggest that phase shift duration is directly 
proportional to the temporal compactness or density of activated 
synapses, because, bursting of in-phase action potentials results 
in an average synaptic potential shift in the frequency of neurons 
and consequently phase shift duration is expected to be inversely 
related to neural density, that is, the higher the packing den- 
sity then the shorter the phase shift duration. This is like local 
neighbors quickly communicating whereas long distance cousins 
take more time and effort to synchronize (i.e., long phase shift 
durations). 

A series of EEG PR studies by Lehmann et al. (2006) and 
Thatcher et al. (2008a, 2009a,b) are consistent with physiologi- 
cal models of EEG PR and have added to earlier studies by Varela 
(1995), Breakspear (2002, 2004), Freeman (2003), and Freeman 
et al. (2003) by measuring discontinuities of electrical potentials 
and current sources of the two main physiological processes that 
underlie PR, namely, phase shift followed by phase lock. Lehmann 
et al. (2006) and Thatcher et al. (2007) demonstrated tempo- 
ral discontinuities of EEG current sources of about 40-250 ms. 
Thatcher et al. (2009a) studied the development of scalp electrode 
distance and PR times by measuring phase shift durations (range 
of about 30-70 ms) and phase lock durations (100-800 ms) from 
birth to 16 years of age where short distance inter-electrode pair- 
ing (6 cm) exhibited shorter phase shift duration and longer 
phase lock duration than longer distance inter- electrode parings 
(18-24 cm). Furthermore, it has been shown that phase shift 
duration is positively related to intelligence while phase lock dura- 
tion is negatively related to intelligence measured by WISC-R 
I.Q. test (Thatcher et al., 2008a). The findings of Thatcher et al. 
(2008a, 2009a) and Lehmann et al. (2006) are consistent with the 
hypothesis that phase shift is a process involved in the recruit- 
ment of available neurons at a given moment of time and phase 
lock duration is the binding or synchrony of groups of neu- 
rons that simultaneously mediate different functions in different 



brain regions (i.e., sustained commitment of neurons). It is also 
possible that phase lock reflects the inhibition of billions of "irrel- 
evant" neurons that are excluded or restricted resulting in the 
"protection" of a small subset of neural loops that mediate a mul- 
tidimensional sub-network. That is, the large spatial inhibition 
isolates a spatially smaller subset of synchronized neurons that 
are masked or invisible to the scalp recorded EEG. The present 
study is a further exploration of the scalp surface EEG studies of 
phase shift and phase lock duration by applying 4-dimensional 
neuroimaging (tEEG) of current sources using Low Resolution 
Electromagnetic Tomography (LORETA) using the time series of 
the center voxel of each of 88 Brodmann areas that comprise 
the DMN (Pascual-Marqui et al, 1994; Pascual-Marqui, 1999; 
Lehmann et al., 2006; Canuet et al., 2011; Langer et al, 2011). 
The present study is designed to explore the nature of sudden 
phase shifts followed by phase lock in small 3-dimensional vol- 
umes of EEG current density located in the center of Brodmann 
areas that constitute the DMN (Buckner et al., 2008). Because of 
the large number of possible network combinations and frequen- 
cies, we limited this study to the delta frequency band (1-4 Hz) 
and the DMN. Analyses of different frequency bands and loca- 
tions show similar basic time domain measures with evidence of 
spatial-frequency "preferences." These analyses will be published 
in the future. 

METHODS 
SUBJECTS 

A total of 70 subjects ranging in age from 13.01 to 19.98 years 
(males = 41) were included in this study. The subjects in the 
study were recruited using newspaper advertisements in rural and 
urban Maryland (Thatcher et al, 1987, 2003, 2007). The inclu- 
sion/exclusion criteria were no history of neurological disorders 
such as epilepsy, head injuries and reported normal develop- 
ment and successful school performance. None of the subjects 
had taken medication of any kind at least 24 h before testing. 
All of the subjects were within the normal range of intelligence 
as measured by the WISC-R and were performing at grade level 
in reading, spelling and arithmetic as measured by the WRAT 
and none were classified as learning disabled nor were any of 
the school aged children in special education classes. All subjects 
were given an eight-item "laterality" test consisting of three tasks 
to determine eye dominance, two tasks to determine foot dom- 
inance, and three tasks to determine hand dominance. Scores 
ranged from —8 (representing strong sinistral preference or left 
handedness), to +8 (representing strong dextral preference or 
right handedness). Dextral dominant children were defined as 
having a laterality score of >2 and sinistral dominant children 
were defined as having a laterality score of < —2. Only 9% of the 
subjects had laterality scores < —2 and 87% of the subjects had 
laterality scores >2 and thus the majority of subjects were right 
side dominant. 

EEG RECORDING 

The EEG was recorded from 19 scalp locations based on the 
International 10/20 system of electrode placement, using linked 
ears as a reference. Eye movement electrodes were applied to 
the inner and outer canthus to monitor artifact and all EEG 
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records were visually inspected and manually edited to remove 
any visible artifact. Two 5 min of EEG was recorded in the eyes 
closed and in the eyes open condition. The order of record- 
ing for the eyes open followed by closed conditions and vice 
versa was counter-balanced across subjects. Each EEG record was 
plotted and visually examined and split-half reliability and test- 
retest reliability measures of the artifacted data were computed 
using the Neuroguide software program (NeuroGuide, v2.6.9). 
The amplifier bandwidths were nominally 1.0-30 Hz, the outputs 
being 3 db down at these frequencies and the EEG was digitized 
at 100 Hz. Analyses were performed on 58 s to 2 min 17 s seg- 
ments of EEG. Split-half reliability tests were conducted on the 
edited EEG segments and only records with > 90% reliability were 
entered into the spectral analyses. Phase shift and lock duration 
were computed only on contiguous EEG segments. 

L0RETA TIME DOMAIN COMPUTATION 

The standard procedures for the computation of LORETA 
were followed according to Pascual-Marqui et al. (1994, 2001), 
Pascual-Marqui (1999), Gomez and Thatcher (2001), Thatcher 
et al. (2005a,b). Numerous studies have demonstrated that 19 
scalp electrodes are sufficient in number to measure intracra- 
nial sources including from the hippocampus (see the listing 
of 795 publications at: http://www.uzh.ch/keyinst/NewLORETA/ 
QuoteLORETA/PapersThatQuoteLORETA05.htm). 

The Talairach Atlas coordinates of the Montreal Neurological 
Institute's MRI average of 305 brains (Pascual-Marqui, 1999; 
Lancaster et al., 2000) was used and the linkage to the standard 
anatomical 7x7x7 mm voxels each with a distinct Talairach 
Atlas Coordinate. Groups of voxels are also defined by the clear 
anatomical landmarks established by Brodmann in 1909 and 
referred to as Brodmann areas. The time series of current source 
vectors in the x, y, and z directions were computed at the center 
voxel of each of 14 Brodmann area that comprises the DMN. In 
addition to the x, y, and z time series from each voxel the resultant 
vector was computed as the square root of the sum of the squares 
for the x, y, and z source moments. 

HUBERT TRANSFORM AND COMPLEX DEMODULATION 

The Hilbert transform of the LORETA time series was computed 
using complex demodulation to compute instantaneous coher- 
ence and phase-differences between each pair of the Brodmann 
area time series with Talaraich atlas coordinates described in 
Table 1 (Granger and Hatanka, 1964; Otnes and Enochson, 1978; 
Bloomfield, 2000). A total of 91 pairs of the LORETA Brodmann 
area time series were used to compute "instantaneous" phase dif- 
ferences. This method is an analytic linear shift-invariant trans- 
form that first multiplies a time series by the complex function 
of a sine and cosine at a specific center frequency (Center fre- 
quency = 2.5 Hz) followed by a low pass filter (6th order low-pass 
Butterworth, bandwidth = 1-4 Hz) which removes all but very 
low frequencies (shifts frequency to 0) and transforms the time 
series into instantaneous amplitude and phase and an "instanta- 
neous" spectrum (Bloomfield, 2000). We place quotations around 
the term "instantaneous" to emphasize that there is always a 
trade-off between time resolution and frequency resolution. The 
broader the band width the higher the time resolution but the 



Table 1 | List of the Brodmann areas of the Default Mode Network 
(DMN) based on Buckner et al. (2008). 
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The Key Institute's Talairach atlas coordinates in the x, y and z directions from 
Lancaster et al. (2000) and Pascual-Marqui (2004). These coordinate values were 
used to calculate the Euclidean Distances as applied to Equation (1) and in 
Figure 5 



lower the frequency resolution and vice versa. Mathematical 
details are in Thatcher et al. (2008a). 

COMPUTATION OF THE 1ST AND 2ND DERIVATIVES OF THE TIME 
SERIES OF PHASE DIFFERENCES 

The 1st and 2nd derivatives of the time series of instantaneous 
phase-differences was calculated for all pair wise combinations of 
DMN voxels in the x, y, and z direction in order to detect instan- 
taneous advancements and reductions of phase-differences. The 
same mathematical procedures published for measuring phase 
shift and phase lock duration of the scalp surface EEG were used 
for the computation of PR of the LORETA time series (Thatcher 
et al, 2008a, 2009a,b). 

PR is composed of two events: ( 1 ) a phase shift of a finite dura- 
tion (SD) and (2) followed by an extended period of phase locking 
as measured by the phase lock duration (LD) and PR = SD + LD. 
Phase Shift duration (SD) is the interval of time from the onset 
of a phase shift to the termination of phase shift (5° thresh- 
old) where the termination is defined by two conditions: (1) a 
peak in the 1st derivative (i.e., 1st derivative changes sign from 
positive to zero to negative) and (2) a peak in the 2nd deriva- 
tive or inflection on the declining side of the time series of first 
derivatives. The peak of the 2nd derivative marked the end of 
the phase shift period. Phase shift duration is the difference in 
time between phase shift onset and phase shift offset or SD(t) = 
S(f)onset — offset- Phase lock duration (LD) was defined as the 
interval of time between the end of a significant phase shift (i.e., 
peak of the 2nd derivative) and the beginning of a subsequent 
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FIGURE 1 | Example of phase reset of LORETA current density time 
series from one subject. Top are the LORETA EEG phase differences with 
respect to the left Hemisphere Brodmann area (BA) 8 time series. The last 
four traces are the phase difference (degrees) for BA8L-BA28L, BA8-BA35L, 



BA8L-BA36L, and BA8L-BA37L. Bottom are the 1st derivatives of the phase 
differences in the top traces in degrees/centiseconds. A 1st derivative 
> 5°/cs marked the onset of a phase shift and an interval of time following 
the phase shift where the 1st derivative ~0 defined the phase lock duration. 



significant phase shift, i.e., marked by the peak of the 2nd deriva- 
tive and the presence of a peak in the 1st derivative or LD(f) = 
S(f)offset — 5(f) onse t- In summary, two measures of phase dynam- 
ics were computed: (1) Phase shift duration (ms) (SD) and, (2) 
Phase lock duration (ms) (LD). Given the range of epoch sizes 
from 58 s to 2 min the range of PRs per subject was 66 to 457. 
Figure 1 shows an example of the computation of PR metrics in a 
single subject. 

DEFAULT MODE NETWORK (DMN) 

Complex Demodulation was used to compute the Hilbert trans- 
form of the current source density time series from 14 left and 
14 right hemisphere Brodmann areas were selected based on 
the review of the human DMN by Buckner et al. (2008). The 
mathematical description of the equivalence of complex demodu- 
lation and the Hilbert transform is by Bloomfield (2000). Table 1 
shows a listing of the nearest Brodmann area fit to the center 
of the Brodmann areas comprising the DMN to the LORETA 
Talairach Atlas coordinates (Pascual-Marqui, 1999; Lancaster 
et al., 2000) and the linkage to the standard anatomical 7 x 
7x7 mm voxels in the approximate center of each Brodmann 
areas. All pair wise combinations of the 14 DMN Brodmann areas 
produced 91 pairs from the left and 91 pairs from the right hemi- 
sphere. Only intra-hemisphere analyses of the delta frequency 
band (1-4 Hz) were included in this study. Different frequen- 
cies and cross-hemisphere combinations will be analyzed in a 
future study. Table 1 shows the Brodmann areas and Talairach 



Atlas coordinates of the voxels in the present study which are 
used to calculate distances between Brodmann area center vox- 
els including the computation of the Euclidean distance between 
Brodmann areas: D = ^Jx 2 + y 2 + z 2 . 

Figure 2 shows the LORETA saggital, coronal and horizontal 
sections of the DMN voxels used in this study. 

RESULTS 

VARIMAX FACTOR ANALYSIS 

A principal components analysis followed by a varimax rotation 
were computed for the 91 left and 91 right intra-hemisphere 
Brodmann area combinations that comprised the DMN for both 
phase shift and phase lock duration. A varimax rotation was used 
because it uses a min/max method to maximize loadings on a 
given component. Table 2 shows that with an eigenvalue cutoff 
of 1.0 that phase shift involved more factors than did phase lock 
duration and that the same number of factors were involved for 
the left and right hemispheres. The results of the factor analysis 
also demonstrated well ordered orthogonality or independent 
clusters of phase shift and phase lock duration where the variable 
loadings were essentially the same in the x, y, z directions and for 
the resultant vector within each hemisphere. 

DISCRETE TEMPORAL FRAMES OR "QUANTA" OF PHASE SHIFT 
DURATION 

Examination of the variables in the factor analyses showed that 
there are three main "modes" or time frames of duration for 
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FIGURE 2 | The locations of the six Default Mode Networks as summarized by Buckner et al. (2008) and represented by the Key Institute LORETA 
voxels (Lancaster et al., 2000; Pascual-Marqui, 2004). 



Table 2 | Number of factors that had Eigenvalues of 1.0 or greater for 
LORETA phase reset in 91 combinations of Brodmann areas for phase 
shift and lock in the 3-dimensions (x, y, z) and also the resultant 
vector R = <J x 2 + y 2 + z 2 in the left and right hemispheres. 

Number of factors phase shift and phase lock 
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Phase lock 





Left 


Right 


Left 


Right 


X 


18 


20 


11 


11 
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19 
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19 
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18 


18 


12 


13 
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18.5 


18.5 


11.75 


12.5 



phase shift duration. Figure 3 shows phase shift duration in 
the x-axis and the percentage of subjects in the study showing 
specific phase shift durations on the y-axis. The x, y, and z direc- 
tions all showed essentially the same phase shift duration modes. 
The three duration modes or "quanta" for phase shift using the 
resultant vector for both the eyes closed and eyes open condi- 
tions demonstrated that discrete durations are present and that 
there is no or minimal overlap between modes. It can be seen in 



Figure 3 that the subjects clustered in three different and discrete 
time frames or modes. Mode 1 showed a peak phase shift dura- 
tion of approximately 25 ms, Mode 2 was approximately 50 ms 
and Mode 3 was approximately 65 ms. There was little overlap 
between phase duration modes and all variables exhibited only 
one of the three modes which demonstrates discrete time frames 
of phase shift duration in specific groups of Brodmann areas that 
comprise the DMN. 

DISCRETE TEMPORAL FRAMES OR "QUANTA" OF PHASE LOCK 
DURATION 

Examination of the factor analyses showed that there are two 
main "modes" or time frames of duration for phase lock duration. 
The x, y, and z directions all showed essentially the same duration 
modes and Figure 4 shows the two duration frames for phase lock 
using the resultant vector for both the eyes closed and eyes open 
conditions. The x-axis is phase lock duration in milliseconds and 
the y-axis are the percent of subjects (N = 70) in both the eyes 
closed and open conditions. It can be seen in Figure 4 that the 
subjects clustered in two different and discrete time frames or 
modes. Mode 1 showed a peak phase lock duration of approx- 
imately 250 ms and Mode 2 was approximately 425 ms. There 
was a minor mode at approximately 800 ms which was much 
smaller than modes 1 and 2. There was little overlap between 
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FIGURE 3 | Phase shift durations between Brodmann areas in the x, y, 2 
LORETA time series directions and the resultant vector in the lower right 
where R = -J x 2 + y+z 2 . The x-axis is phase shift duration in milliseconds 
and the y-axis is the percent of subjects that exhibited a given phase shift 
duration for different Brodmann area pairs. The solid line is the eyes closed 
condition and the dashed line is the eyes open condition. All of the subjects 



are represented within each curve. For example, 100% of the subjects 
exhibited a phase shift duration between 1 8 and 35 ms for Brodmann areas 8 
and 9 (upper left panel x-shift) and similarly for each Brodmann area pair. The 
finding of discrete phase shift durations with none or little overlap of data 
points under each phase shift duration curve was a dominant feature of 
phase shift duration and demonstrates discrete "temporal quanta." 



phase duration modes and all variables exhibited only one of the 
two modes with no bi-modal distributions which demonstrates 
discrete time frames of phase lock duration in specific groups of 
Brodmann areas that comprise the DMN. 

SHORT vs. LONG DISTANCE CONNECTIONS AND LORETA PHASE RESET 

The finding of discrete time frames or modes was explored fur- 
ther by correlating the Euclidean distance of the Talarich atlas 
coordinates of the DMN (i.e., square root of the sum of squares 
of x + y + z) with phase shift and lock duration. Figure 5 are 
the results of regression analyses in which statistically significant 
inverse relationships were present between phase shift vs. phase 
lock duration. That is, phase shift duration tended to be short 
when Brodmann areas were closer together and lengthened as the 
distance between Brodmann areas increased. In contrast, phase 
lock duration tended to be long when Brodmann areas were closer 
together and short when Brodmann areas were more distant. 

These analyses showed that both a non-linear relationship 
involving "modes" or discrete time frames were present if a small 
number of Brodmann areas are compared as well as a relatively 
smooth exponential decline as a function of distance between 
Brodmann areas when all 91 Brodmann areas are examined. 



In order to further explore both the general and discrete aspect 
of the findings we evaluated each combination of Brodmann areas 
by fitting an exponential equation to the data points. Exponentials 
provided the best fit of the data, for example, the red line is the 
evaluation of the equation: 

T=b l +e V / (1) 

where T = duration (ms), d = distance between Brodmann areas 
(millimeters) and b\, b^, and bj, are coefficients. The coefficients 
were derived using Data Desk statistical package (Velleman, 
1997). The regression to Equation (1) was statistically significant 
but in opposite directions for phase shift vs. phase lock duration. 
Brodmann area pairs that exhibited short phase shift durations 
also exhibited long phase lock durations while Brodmann areas 
that exhibited long phase shift durations exhibited short phase 
lock durations as shown in Figure 5. The coefficients of the equa- 
tion are in Table 3. The coefficients in Table 3 allow one to test the 
hypotheses in this study by simply entering Talairach Euclidian 
distances between voxels of the brain into an Excel worksheet and 
then compute the predicted phase lock and phase shift durations. 
Shifts in the intercept b\ will result in shifts in the starting point 
of the function. 
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FIGURE 4 I Phase lock durations between Brodmann areas in the x, y, z 
LORETA time series directions and the resultant vector in the lower right 
where R = -J x 2 + y+z 2 . The x-axis is phase lock duration in milliseconds and 
the y-axis is the percent of subjects that exhibited a given phase lock duration 
for different Brodmann area pairs. The solid line is the eyes closed condition 
and the dashed line is the eyes open condition. All of the subjects are 



represented within each curve. For example, 100% of the subjects exhibited 
a phase shift duration between 250 and 500 ms for Brodmann areas 8 and 10 
(upper left panel x-lock) and similarly for each Brodmann area pair. The finding 
of discrete phase lock durations with none or little overlap of data points 
under each phase lock duration curve was a dominant feature of phase lock 
duration and demonstrates discrete "temporal quanta." 



These data show that nearest neighbors exhibit short phase 
shift and long phase lock durations while long distance rela- 
tions exhibit long phase shift and short phase lock durations. 
This is analogous to every day life where one quickly engages 
and maintains communication with a local neighbor but it 
takes longer to engage and remain connected to a long distance 
neighbor. 

EYES OPEN vs. EYES CLOSED CONDITIONS AND LEFT vs. RIGHT 
HEMISPHERE 

The eyes open condition was an independent replication of the 
eyes closed condition because the eye conditions were recorded at 
different times and the order of recording was counter balanced 
across subjects. Multivariate analysis of variance (MAN OVA) with 
Bonferroni correction was performed where eyes open and eyes 
closed conditions and left vs. right hemisphere as factors for all 91 
Brodmann area combinations. The MANOVA was used because 
it controls for the intercorrelation between the dependent and 
between the independent variables. Table 4 shows the results of 
these analyses which showed that eyes open vs. eyes closed con- 
ditions were not statistically significantly different in the x, y, z 



directions and resultant vector for phase shift duration. However, 
eyes open vs. closed conditions were statistically significant for 
phase lock duration in the x and z directions and for the resul- 
tant vector. In all cases, phase lock duration in the eyes closed 
condition was longer than in the eyes open condition (mean 
difference = 9.52 ms). 

As seen in Table 4, there were no statistically significant differ- 
ences between left vs. right hemisphere in phase shift duration, 
however, there were statistically significant differences in phase 
lock duration in the x and y directions. In all cases, phase lock 
duration was longer in the right hemisphere than in the left 
hemisphere (mean = 11.05 ms). 

VOLUME CONDUCTION TEST 

Finally, the role of volume conduction is important to evalu- 
ate. Volume conduction is the propagation of an electomagnetic 
field at the speed of light or 3 x 10 10 cm/s (Feynman et al., 
1963; Malmivuo and Plonsey, 1995). For the distance of the 
order of centimeters the delay is approximately 3.3 x s 
which is a delay that is extremely small and therefore approx- 
imates zero phase difference at all points in a volume. Phase 
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FIGURE 5 | The x-axis is the Euclidean distance between the center 
voxels that comprise the DMN Brodmann areas as described in Table 1. 

The y-axis is phase shift duration (Top) and phase lock duration (Bottom) of 
the Brodmann areas described in Table 1. The left row are the left 
hemisphere Brodmann areas and the right row are the right hemisphere 
Brodmann areas. The red line is the fit of an exponential equation 



T = fr| + gte+fa/d) where T = duration time (ms), d = distance between 
Brodmann areas (mm) and h\, £>2. and £13 are coefficients. Ft = regression 
correlation and p = statistical probability. Phase shift and phase lock are 
inversely related where Brodmann areas with short phase shift duration 
exhibit long phase lock durations while Brodmann areas with short phase lock 
durations exhibit long phase shift durations. 



difference is measured by the "imaginary number" component of 
the cross-spectrum and PR is the 1st derivative of the "imaginary 
number" that by definition is not volume conduction if a phase 
difference is greater than zero. A test of the possible influence of 
volume conduction involved measuring the absolute phase differ- 
ence between Brodmann areas and then determining if the phase 
difference was greater than zero and how phase difference changes 
with distance. If greater than zero then the phase lock measures 
cannot be explained by volume conduction. Further, because 
connection density decreases as a function of distance and con- 
duction delays linearly accumulate with distance then if there is 
an increase in phase differences as a function of distance than 
this also cannot be explained by volume conduction. If absolute 
phase differences are significantly greater than zero and increase 



as a function of distance than this is consistent with physiological 
connections between connected parts of a network in which con- 
duction velocity, synaptic rise times and synaptic delays produce 
accumulative time delays. 

Figure 6 shows the results of this test in which mean phase dif- 
ferences were not equal to zero and instead varied as a function 
of distance which is what is expected in a network of connections 
and cannot be explained by volume conduction. In fact, phase 
difference values were many times greater than zero and even the 
most close or nearby neighbor Brodmann areas 8 and 9 differed 
by 2-4° and therefore also cannot be explained solely by vol- 
ume conduction. Further, phase lock duration is maximal at short 
distances and declines with distance that also cannot be explained 
by volume conduction because a sustained phase difference over 
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time must be maintained which is consistent with the known 
physiology of local and distant loop iteration that can extend for 
many milliseconds (Beste and Dinse, 2013). 

DISCUSSION 

This study extends the investigation of the spatial and temporal 
properties of the scalp surface EEG PR to the current sources 
of the EEG. There was an approximate seven year age range of 

Table 3 | Coefficients of equation (1) {bi, bz, and 03) for the resultant 
vector for phase lock and phase shift durations in the left and right 
hemispheres used in Figure 5. 



b1 


b2 


b3 




PHASE SHIFT DURATION 


Left -65.0023 
Right -49.7859 


4.9086 
4.7859 


-3.2892 
-3.7078 




PHASE LOCK DURATION 


Left -781.9091 
Right -1521.5681 


6.8430 
7.4068 


5.9441 
4.5508 


Table 4 | The results of the MANOVA for eyes closed vs. eyes open 
conditions and hemisphere for phase lock and phase shift durations. 


Phase SHIFT 


Phase LOCK 


F-ratio P-value 


F-ratio 


P-value 



EC vs. EO 

X 0.2335 
Y 1.3148 
Z 0.0623 
R 0.0123 



ns 
ns 
ns 
ns 



11.8840 
2.7337 
3.8381 
9.8722 



0.0006 

ns 
0.0501 
0.0017 





LEFT vs. RIGHT 


X 0.4029 


ns 


24.0740 


0.0001 


Y 0.3493 


ns 


8.2249 


0.0041 


Z 0.7527 


ns 


1.5469 


ns 


R 2.9703 


ns 


2.9038 


ns 



subjects in this study, however, no significant correlations with 
age were present. A new finding is discrete temporal modes of 
phase shift and phase lock duration that are unique for different 
Brodmann area associations of the DMN. A second finding is sig- 
nificantly greater complexity and higher dimensionality of phase 
shift duration in comparison to phase lock duration, i.e., a tempo- 
ral dimension reduction from phase shift to phase lock duration. 
A third finding is an inverse relationship between spatial distance 
between Brodmann areas and PR metrics where short distances 
are exponentially related to long phase lock durations and short 
phase shift duration, whereas, long distances are exponentially 
related to long phase shift durations and short phase lock dura- 
tion. Overall, the findings are consistent with previous studies 
demonstrating that temporal "frames" or "chunks" are discrete 
and different for different Brodmann areas of the human brain. 
The results can be explained by a local vs. distant connection 
model which is unique for different Brodmann areas and thereby 
accounting for discrete temporal quanta while also exhibiting an 
exponential decrease in the number of connections with distance 
that corresponds to exponential changes in PR as a function of 
distance. Furthermore, the data indicate that discrete peaks in 
PR duration reflect a different or unique ratio of long vs. short 
distance connection densities within different Brodmann areas. 

Studies by Ko and Ermentrout (2007) and Tiesinga and 
Sejnowski (2010) use mathematical models that best explain the 
onset of a phase shift to be initiated by cortico-cortical long dis- 
tant excitatory dendritic synapses. The mathematical models and 
experiments fit the physiological facts that local distance con- 
nections are exponentially more numerous and exhibit temporal 
compactness while distant connections are less numerous and 
exhibit temporal dispersion. Therefore, the ratio of local to dis- 
tant connections determines the duration of phase shift and this 
ratio varies as an exponential function of distance between a given 
Brodmann area. This conclusion is consistent with studies show- 
ing anti-phase shifts are related to neural packing density (Li 
and Zhou, 2011) and with the surface EEG PR where increased 
local packing density was hypothesized to explain the difference 
between PR in the posterior- to-anterior direction (e.g., 01-P3) 
vs. the anterior-to-posterior directions (e.g., Fpl-F3) because of 



Correlations of Default Brain Phase Difference with Euclidean Difference Distances 
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FIGURE 6 I Test of volume conduction. The y-axis is the mean absolute phase differences (degrees) of the resultant vector between Brodmann areas. The 
x-axis is the Euclidian distance (mm) between all Brodmann area pairs. The left graph is from the left hemisphere and the right graph is from the right hemisphere. 
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increased packing density in the posterior cortex as compared to 
frontal cortex (Thatcher et al., 2009a). The results of the present 
study are also consistent with the hypothesis that phase shift is 
related to recruitment of neurons and phase lock duration is 
directly related to synchrony or binding of neurons as an expo- 
nential function of packing density and inter-node (Brodmann 
area) distance. 

LIMITATIONS OF THIS STUDY 

One limitation is the sample rate used where 100 Hz sample 
rates limit analyses to 10 ms resolution and therefore phase shift 
durations at shorter durations may be present but were missed. 
However, the time resolutions appears to be due more to phys- 
iology than the speed of computer measurement because we 
have measured phase shift and lock duration at sample rates 
of 512 Hz (1.95 ms resolution) and find the same minimum 
phase shift durations of about 20 ms as reported in this study. 
Mathematically the limiting relationship is between distance and 
time as described by Equation (1) that can be evaluated using 
the coefficients in Table 3. If d = 0 then duration = infinity. If 
d = infinity then duration = 0. The values between zero and 
infinity are a good fit of how packing density of neurons as a 
function of distance are related to PR duration. Higher sam- 
ple rates should be used to further test the temporal synchrony 
limits as reported in this study. Another limitation are analyses 
only of the delta frequency band. Because of the large volume 
of data and limitations of publication space it was necessary to 
limit the analyses to a single frequency band, i.e., the delta fre- 
quency band (1-4 Hz). Different relations are present at different 
frequency bands as well as cross-frequency band coupling and 
these relations have been measured and will be the subject of 
future publications. Another limitation is the use of a center 
frequency of a narrow band that is necessary with the Hilbert 
transform. The Gabor transform is independent of a center fre- 
quency and provides optimal time frequency resolution (Witte 
and Schack, 2003) and higher temporal-spatial resolution and can 
yield more detailed time-frequency information than Complex 
Demodulation. While cross-frequency phase shift and phase lock 
measurements is important, however, because of page limitations 
this topic will be also presented in future studies. 

DISCRETE DURATIONS AND THE DEFAULT MODE NETWORK 

Shallace (1964), Efron (1967, 1970a,b), Allport (1968), and oth- 
ers (Sanford, 1971; Varela, 1995; Varela et al., 2001) have shown a 
minimum perceptual frame from approximately 40 ms for audi- 
tory stimuli to 140 ms for visual stimuli that are durations that 
temporally distinguish events as being successive in time where 
duration is defined at Ti - T2 = 0, i.e., simultaneity where there 
is no perceived time difference between two distinct events. These 
studies as well as others show that learning-dependent changes in 
neural networks is not a continuous process but rather a discon- 
tinuous sequencing of narrow time windows (Thatcher and John, 
1977; John, 2005; Lehmann et al, 2006; Thatcher et al, 2007, 
2008b, 2009a). Thatcher et al. (2009a) indicated a linkage between 
spontaneous and ongoing perceptual frames and event related 
desynchronization (ERD) by considering phase shift duration and 



phase lock duration as elemental "atoms" that underlie the dura- 
tion of perceptual frames and ERD. For example, in Thatcher 
et al. (2009a) the mode of scalp surface EEG PR was tempo- 
rally bounded with a minimal phase shift duration of about 45 ms 
and a maximum phase shift duration of about 70 ms. Phase shift 
was followed by phase lock that was temporally bounded from 
about 150 ms to about 800 ms with the most frequent phase lock- 
ing intervals between 200 and 450 ms (Thatcher et al., 2009a). 
The findings in the present study are consistent with the ear- 
lier surface EEG analyses and indicate that 3-dimensional current 
source phase shift and phase lock are ongoing spontaneous pro- 
cesses that occur between network nodes and that discrete phase 
shift durations (i.e., discrete quanta of time) operate at higher 
temporal-spatial precision than the surface EEG. Network nodes 
are defined as clusters of neurons connected to other clusters 
(nodes) and the present study shows that a fundamental prop- 
erty of nodes or clusters is to operate like temporal "shutters" 
that open and close at specific durations. Each Brodmann area 
maintains a different ratio of local vs. distant connections but 
nonetheless follows the general rule of an exponential decrease in 
local connections as a function of distance. This anatomical fact 
may explain the findings of discrete phase shift and lock duration 
due to unique local vs. distant excitatory connection densities in 
different Brodmann areas (Sporns, 2011). 

Figure 7 illustrates an hypothesis to explain the findings in this 
study by fitting the data to a single exponential model based on 
the ratio of local and distant excitatory dendritic synapses. Table 3 
provides the coefficients of Equation ( 1 ) to allow one to experi- 
ment with different Talariach distances, e.g., enter the x, y, and z 
Euclidean distance between Hagmann et al. (2008) Modules in 
millimeters into Equation (1) and then calculate the predicted 
LORETA phase shift and phase lock duration. 

The EEG is the summation of LFPs therefore phase shifts are 
necessarily related to connection density. However, the arrival of 
distant synaptic action potentials exhibits temporal dispersion 
whereas local excitatory connections are higher in number and 
temporally compact. Therefore, it is hypothesized that the ratio of 
distant to local connections varies as a function of distance from 
any Brodmann area and as a consequence longer phase shift dura- 
tion occurs as the ratio shifts toward distant excitatory inputs. It 
is hypothesized that phase lock is inversely proportional to phase 
shift duration based on the spatial-temporal GABA connections 
and delays between excitatory neurotransmitter EPSPs (and exci- 
tatory neuromodulators) in local and long distance loops. It is 
also hypothesized that the physiological differences in the gene- 
sis of phase shift vs. phase lock is related to the "Gap" of time 
between distant and local EPSP excitation that produce the phase 
shift followed by long duration inhibitory synaptic potentials that 
contribute to phase lock duration. The collective resonance of 
the rebound from local inhibition and the arrival of long distant 
EPSPs contribute to the onset of the PR process. 

Whether or not the phase shift and phase lock are time 
correlated to a task is irrelevant since "self-organized critical- 
ity" is an ongoing background emergent process that on the 
average produces an approximately 20-80 ms period of phase 
shift or "uncertainty" or approximate duration of "chaos" fol- 
lowed on the average by a 200-800 ms period of phase locking 
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FIGURE 7 | (A) (Top) is an example of discrete durations or "tempora 
quanta" of phase shift duration in different Brodmann areas of the DMN. The 
eyes vs. eyes closed shift to longer durations are due to increased functional 
connectivity with increased input that is like a "shutter" whose duration is 
proportional to the number of recruited neurons. The discontinuities are due 
to different packing densities in different Brodmann areas. The bottom plot 
(B) is the evaluation of Equation (1) using the mean LORETA phase shift and 
lock duration for both eyes closed and open conditions for the 91 Brodmann 
area combinations: T = fr| + e b2+ ( b3/d 'l The link between EEG phase shift 
and neural packing density is the physiological observation of action potential 
bursts when in-phase to LFPs vs. suppression of action potentials when 
ant-phase to LFPs Hughes and Crunelli (2007). EEG is the summation of 
LFPs therefore phase shifts are necessarily related to neural packing density, 
i.e., the higher the packing density than the longer the phase shift as a 
property of summation. It is hypothesized that the increased shift duration 



between eyes closed vs. eyes open is due to increased arousal and 
increased depolarization resulting in increased functional connectivity 
(increased neural resource). Phase lock is inversely proportional to phase shift 
duration based on the spatial-temporal GABA connections and delays 
between excitatory neurotransmitter EPSPs in local and long distance loops. 
The physiological differences in the genesis of phase shift vs. phase lock is 
related to the "Gap" of time which is a transition time between local EPSP 
excitation and re-enterant long distance EPSPs that produce the phase shift 
followed by long duration inhibitory synaptic potentials that contribute to 
phaselock duration. The rebound from inhibition and arrival of EPSPs starts 
the phase reset process. Predicted phase shift and lock durations can be 
evaluated by using the coefficients in Table 3. For example, the distance 
between Modulel and Module2 from Hagmann et al. (2008) and Thatcher 
et al. (2011 , Table 3) is 43.4 mm which that predicts a phase shift duration = 
56 ms and phase lock duration = 300 ms based on Equation (1). 



or "stability." This process represent a continuous sequence of 
meta-stable states (Rabinovich et al., 2012). The loop network 
background process includes "blank" periods when large assem- 
blies of neurons are in a PR mode (i.e., phase shift and phase 
locking) and otherwise not available to participate in loops which 



represents an approximate average 135 ms "Gap" or period of 
time between when neurons recruited by a phase shift are fol- 
lowed by phase lock onset. This gap is a statistical region with an 
average and range between the end of a phase shift and the begin- 
ning of a subsequent phase lock. This appears to be a fundamental 
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"unconscious" transition time that is near to the flicker fusion fre- 
quency and may be why TV viewers rely upon instant replay to 
confirm a touch down in football or a replay in baseball, etc. (i.e., 
human time frames are too long and require higher speed instant 
replay to determine the reality of an event under question). The 
"gap" interval of about 135 ms is spatially distributed and tempo- 
rally averaged across billions of loops in the brain and results in a 
continuous or smooth appearance of reality. 

Another type of "refractory" period or "gap" is when phase 
locked neurons are unavailable for allocation by a different clus- 
ter of neurons at a different moment of time. Long distance phase 
locking of local clusters of neurons can result in a reduction in 
the amplitude of the surface EEG because phase locking occurs 
over long distances and thus reduces the size of the number 
of synchronized local cluster of neurons by spatial differentia- 
tion. This is consistent with studies of schizophrenia that show 
hyperconnectivity in local frontal and parietal regions associated 
with increased local current density (Canuet et al., 2011). The 
hypothesis of linking PR during the background spontaneous 
EEG provides a new definition of the term "desynchronization" 
used to describe ERD and the waxing and waning of the spon- 
taneous EEG. That is desynchronization is actually "spatially 
differentiated phase reset" or "micro bonding" of local clusters 
of neurons connected across long distances to other local clusters 
for brief periods of time (Thatcher et al., 2009a). 

EYES OPEN vs. CLOSED AND LEFT vs. RIGHT HEMISPHERE 

The eyes closed vs. eyes open conditions were measured at differ- 
ent times with a pause between recordings and thus are within 
subject replications. Both conditions exhibited significant fits 
to the exponential Equation (1) and there were no significant 
differences between eyes open and eyes closed phase shift dura- 
tion (see Table 4). The eyes closed condition did exhibit a small 
but significantly longer phase lock duration which may reflect 
the relatively small differences in neural excitability between these 
two states. The lack of a large difference between resting eyes 
closed and eyes open states is consistent with fMRI studies of 
the DMN that is evident in the resting state (Raichle, 2010). This 
is because the resting state is characterized by quiet repose with 
either eyes closed or eyes open and with or without visual fixation. 
During the resting state subjects typically experience an ongoing 
state of conscious awareness filled with "stimulus-independent 
thoughts" (SITs; Antrobus, 1968) or more commonly, daydream- 
ing or mind wandering (Mason et al., 2007). Internal thoughts, 
self-awareness and rumination commonly dominate the resting 
state as opposed to active task engagement (Sridharan et al., 2009; 
Raichle, 2010). The right hemisphere exhibited longer phase lock 
durations than the left hemisphere (see Table 4) which is con- 
sistent with the findings of lower packing density in the right 
hemisphere in comparison to the left, especially in the planum 
temporale (Gur et al, 1980; Buchell et al, 2004). 
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